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A trajectory in the Schrodinger wave for an electron in an attractive Coulomb potential with 
the dynamical behavior is proposed and illustrated for a scattering and a bound state. The 
scattering cross section derived from the trajectories is almost exactly equal to that from the 
usual wave theory. The statistical nature of the result is examined. The period of one cycle in 
the bound state is exactly equal to that of the corresponding classical motion. 

1 Introduction 

Nowadays, we can move an atom to put on the arbitrary positions on metal. We can measure 
whether one electron adhered to the capacitor. There is no uncertainty to operate one atom. We also 
have ascertained that there is no indeterminacy on determining a submicron diameter by a longer 
wavelength laser precisely. jp It would be one of the urgent requisites to make a dynamical equation 
of motion of an atomic scale 'particle' without any contradiction to the wave theory. 

The trajectory of the harmonic oscillator in the wave equation has been discussed, j^j Since the 
system treats only the bound state, the statistical nature, one of the important characteristics, in 
quantum theory can not be seen clearly. In the present paper, by using the method of the mcf j| |), 
a trajectory of an electron in an attractive Coulomb potential, like the hydrogen atom in a state with 
a positive or negative energy, is proposed. On discussing the cross section for the scattering state, the 
statistical but in principle determinate nature like classical mechanics will be seen as an illustration. 

In section || the dynamical theory that would lead to a trajectory is described. The traveling 
waves and usual stationary wave functions are discussed. In section ^ the trajectory of an electron in 
the scattering state by the Coulomb potential is investigated. In sections [| and || the cross section 
and the flux of the beam of particles are discussed. In section ^ the trajectory of an electron in the 
hydrogen atom in a bound state is analyzed and illustrated. Conclusion and remarks on the wave 
functions are given in section |?|. 

2 Dynamics and wave function 

The dynamics that leads to the mode trajectory of an electron in an attractive Coulomb potential 
with a charge e(> 0) is summarized. || The wave function *S> describing the motion of an electron 
satisfies the Schrodinger equation 

%n y ' = -— A--U(r,t), (1) 



dt \ 2m 

where constant m or — e is electron mass or charge, respectively. 

The equation is assumed to be separable in variables t,xi,X2 and X3. Let the wave function be 

tf(r,t) =e- l£t / R $i(xi,£,a)$ 2 (£2,a,/?)$ 3 (:r3,/3), (2) 
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where E, a and (3 are constants of separation, and E is assumed to be the energy of the system. These 
constants should be called mode parameters. The wave function of the form 



= \ /^r ; iex P {;U' j = 1, 2, 3, (3) 

is sought, where functions Pj's and Wj's are real. This should be called a traveling wave. Let functions 
W/s satisfy the condition that in each classical region of Xj for j = 1, 2, 3 

W j (x j )^W j (x j ) ch (4) 

where the sum of them 

3 

^01 = ^^(^)01 (5) 

J'=l 

is the Hamilton characteristic function of the Hamilton- Jacobi equation in classical mechanics. j| The 
classical region stands for the domain in which the chracteristic function holds true. 
If Wj 's are found uniquely, the sum of them 

W(xi, x 2 , x 3 , E, a,0) = Wx{x x ,E, a) + W 2 {x 2 , a, 0) + W 3 (x 3 ,f3) (6) 

is named the mode characteristic function (mcf) for the system. 
The equations of motion for the electron are assumed 

dW 1, , dW dW 

-dE=n (t - t ^ ^ = c - W = c ^ (7) 

where to, c a and are constants (independent of f) that are determined by initial conditions for the 
system. Variable t is considered the dynamical time for the system. The trajectory derived from 
Eqs. (0) should be called the mode trajectory. 

The square integrable wave function is obtained as follows. Let the endpoints of the Xj coordinate 
be dj and bj. The electron moves in the region dj < Xj < bj. Let it start from a point Xjo to the 
increasing Xj direction and the mcf in the coordinate space be Wj (xj ) . The traveling wave associated 
with the motion reversing from bj to a,j should be assumed to be given by J Pj(xj) exp[i{— Wj(xj) + 
2Wj(bj)}], which is also the solution of Eq. (|l|). 

The wave observed at Xj(> Xjo) should be the superposition of the traveling waves associated with 
the alternating motion of the electron || 

^(xj)-**^ 2 ^. (8) 

This is finite, zero, at bj. If the electron turns at cij and runs to xjq, the mcf should be assumed to 
be given by Wj(xj) — 2Wj(cij) + 2Wj(bj). The wave function for xj(< Xjq) is written as 

- $fcj)e aw 'W + ^ j (x j )e i ^- 2W ^ +2W ^. (9) 

This is finite, zero, at a,j. The wave functions (||) and (^|) are finite for a,j < Xj < bj for any 'mode' 
with real values of parameters. They are not always equal at XjQ. As a wave function, it might be 
multi- valued. If and only if 2[Wj(bj) — Wj(dj)] is a multiple of 2-7T, they are equal and constitute a 
stationay wave function in cij < Xj < bj. 

In what follows, the wave equation ([!]) is analized in the polar coordinate system, (r,6,cp). The 
Hamilton characteristic function is summarized as follows: 

W cl (r, 6, & E, I, ;it) = W r , cl (r, E, I) + W e , c i(9, 1, /*) + pH, (10) 
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where W r>c \ and Wg tC \ are determined from equations 



1 

2m 



'dW, 



r.cl 



V dr 



i 2 h 2 



-- = e, 

r 



dWg,, 

89 



sm 



i 2 h 2 . 



(ii) 

(12) 



Here, E stands for the energy and will be expressed in terms of rj s , (24), for the scattering state or 
r), (|76|), for the bound state. Parameters I and /i are real numbers and lh stands for the angular 
momentum and fih is one of its components. These expressions have been introduced for convenience 
for comparison with the quantum theory. 



3 Scattering state 



The scattering state of an electron in the Coulomb potential is analyzed in the spherical polar coor- 
dinate system. The wave function ^(r, t) is expressed in the spherical polar coordinates with mode 
parameters, E, v and \x as 



*(r,t) = exp(-iEt/?i)^(r,E), 
$(r,P) = R(r,E,v)Y(e,v,fj,)eKp(ifi<j>). 



(13) 
(14) 



Constant E stands for the energy and Tiv for the orbital angular momentum, and hfi represents the 
component of the angular momentum along the polar axis. When v and /i are integral numbers, they 
are usual azimuthal and magnetic quantum numbers as will be seen in Eqs. ( |T5| ) and (|^). p] 

The mcf expressed in terms of the spherical polar coordinates are obtained as follows. The function 
Y(9, i/, /j,) satisfies the differential equation 



d 2 n A r 

+ co t0_ + „(„+l)_ 
dtr dO sin i 



Y(6,v,n)=0. 



(15) 



The solution is a linear combination of linearly independent associated Legendre functions, P^(cos( 
and Q£(cos9). |?J A traveling wave in the 9 coordinate space is given by 



2 . a dW 
-sin0— — 

7T d9 



Q^(cos 

Zr VI 

= yfp e m *. 
The mcf for the 9 component should be determined as 

"ttP£(cos 



-1/2 



exp [iWe(0, u, (j,)] 



We (9, v, /i) = arctan 



2 Q£(cos( 



(16) 



(17) 



because of the similarity to the characteristic function Wg tC \ in the classical region and the validity of 
the results derived from this as will be seen in the following. 

An expression for function tan(Wg) is written for < 9 < tt as [0 



?rP^(cos6 
2Q£(cos6 

OO 

= E 

oo 

E 



(l/2 + M ) fc (l + ;/ + M ) fc . 

— — — - sin U2k + v + [i- 

k\{v + i/2) k 

(1/2 + A») fc (l + f + A*)fc 



fc=0 



fc!(^ + 3/2) fc 



cos[(2fc + i/ + /i + 1) 



(18) 
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The value of We at 9 = or tt is @ 

W 9 (0) = W fl (0,i/,AO=7r/x > W e (7r) = W fl (7r,i/,A*)=7ri/. (19) 

By the asymptotic expansion of the Legendre functions for v 3> 1, @] it can be obtained that 

WM0, i ^)«(^ + y2)0 + ( 1 /4 + 1 /2MK (e < 6 < 7T-e, e > 0). (20) 

Both the derivative cWg/d^ and dWe/dp are monotonic as a function of 9 and very similar to 
those of the characteristic function W$ tC i in a restricted (or classical) region as found by a computer 
calculation. 

Radial wave function satisfies the differential equation 



d 2 v{v+l) 2m /e 2 
dr 2 r 2 fi 2 



where u(r, 2?, v) = rR(r, E, v). 

With E positive the linearly independent solutions are 



u{r) = 0, 



(21) 



UM 
Uy 



= e- lp p v+1 V{v + l + ir ls ,2 V + 2 : i2p) 1 



where 



P = 



2mE 



Vs ~ h \ 2E- 



(22) 
(23) 



(24) 



Function V(a, b, z) is defined for convenience 
V(a,b,z) = r(o) 



U(a, b, z) — cos na — - — M(a, b, z) 



, r(6-o) , , 

= -vrcotTrfe— — -Ma,6,z 

r(&)r(i - a) 
+r(6 - Vjz^Mil + a-b,2-b,z). 



(25) 



Functions M(a, b, z) and t/(a, 6, z) are the Kummer functions. || 

For the far region from the center of the potential, p 3> 1, by putting a = is+l + ii] s and b = 2^ + 2, 
it holds i, Ii0| that 



M (a, b, i2p) 
V(a,b,i2p) 



e tp T{b) 



(2p)- a (l+ 1 j 



r(a*) 



2p 



+ c.c. 



where 



[—i sin(7ra)G(/9, a) — cos(ira)G(p, a) 

ia(l — a*) 



G(p,a) = r(a)e- l(p - 7ra/2) (2 / 9)- a 1 + 



2p 



(26) 
(27) 

(28) 



and G(p,a)* is the complex conjugate (c.c.) of G(p,a). These asymptotic forms indicate that the 
linear combination of functions M and V producing an outgoing traveling wave in the far region from 
the origin should be written as 



u = e- lp p v+1 



V(a, b, i2p) + iM(a, b, %2p) sin(7ra) 



r(o)r(6- a) 



r(6) 



(29) 
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By equations mentioned above, this leads to 

u « e -^ e 1 /2^» 2 - 1 - i '|r(a)| exp[i(p + i ls log(2p) - ar g r(a) - 1/^(1 + i/))] 

^ 1+i Kl±^l_| + 0( ^, (30) 

for p large. J5|, |ll[ It would suggest that the traveling wave in the r coordinate space should be given 
by 

. , .T(a)T(b — a) . , 

u = uy + tUM sin(7nzj — — =e (ur + zuj) 

r(6) 



u R + ujexp I iarctan— ^-J = e^" P r e tWr , (31) 

and thus the mcf in the r coordinate is given by 

W r (r, E, v) = arctan — . (32) 

Ur 

Here, functions ur and uj have been assumed to be real. The parameter v is a real number. The 
ratio ur/uj is written as 

uj cos(tt6) - exp(-27r?7 s ) 

— = ■ , ,s ■«) (33) 

ur sm(7ro) 

Function H is real since it holds by Kummer's transformation 0] that 

e- ip M(a,b,i2p) = e tp M(a,b,i2p)*, (35) 
e~ lp M{a + 1 - 6, 2 - 6, «2p) = e lp Af (a + 1-6,2-6, i2p)* . (36) 

Function (Q with expressions (|l^) and ( |3l| ) specifies a traveling wave associated with the motion 
of an electron in the 'mode' (E, v, p) in the scattering state. 
In the far region from the origin the mcf is approximated as 

W r (r,E,v) ~ p + Vs \og2p - ar g r(a) - Y 27 r(l + v) + " (1 + ^ + ^ 

2p 

+ 0( P - 2 ). (37) 

This is nearly equal to the corresponding Hamilton characteristic function. ||, [ll| It indicates the 
validity of the definition of the mcf (621). For r small, it is obtained that 



W r (0) = W r (0, E, V) = V 2 7T. (38) 

The mcf in the spherical polar coordinates is summarized as 

W(r,e,(f>,E,v,p) = ±W r (r,E,v)±W e (e,p,p)±tMt>, (39) 

where the sign should be adopted according to the direction of the motion of the particle in each 
coordinate. Equations of motion (0) with r, 9, <p, v or p in place of x\, X2, x%, a or /?, respectively, lead 
to a trajectory of the electron. 



■5 



To get the trajectory, it is necessary to calculate derivatives dW r /dE and dW r /dv. It is obtained 
from Eqs. © and (ph that 



dW r 
dv 

d_ ( uj 
dv \Ur 



1 



9 ( ui 



1 + (m//w,r) 2 dv \u R J ' 

exp(-27T77 s ) C0S7T& - 1 
Z7T ^ ±Z 



sin 7r6 



ip(a) - tp(a*) 



2(^(6-l)+^(6)-log2p) 



M 2 [ da db 



1 /<9M SM 
M V~da" + "96" 



where M = M(a,b,i2p),ip(b) is the psi (digamma) function and 

M 2 = M(a + 1 - 6,2- 6,i2p) = M(-i^ + ir) s , -2v,i2p). 



(40) 



(41) 



For example, a trajectory of an electron incident from a point distant from the origin of the 
potential, (po, 9q, 4>q), and scattered to another distant point is considered. To be specific, that 
Po = oo and 9q — is assumed. For the trajectory from po to the origin, the mcf is written as 



W(r, 9, 4>, E, v, p) = -W r +W9+W4,, 
where W$ = pcj). The trajectory is given by the equations (Q), or 



d_ 

dv 



d 



(-W r + W 6 ) 
{We + W+) 



— (-W r (oo) + Wg(0)) = — argr(a)+y 2 7r, 
dv dv 



d 



W e (0) + <h = * + <h- 



dp, dp 
For the path from the origin to p n = p(9 — 7r), the mcf is given by 

W{r, 9, (f>, E, v, p) = W r + W e + W <t> - 2W r (0). 



(42) 

(43) 
(44) 

(45) 



This has been taken to be continuous to the incident mcf, (42), at the origin. 
The trajectory equation is written as 

j- (W r + W e - 2W r (0)) 

= ^argr(a)+y 2 7r, 
dv 



(46) 



and Eq. (|44|). Since function d v W r shows monotonic decrease with respect to p while d u We does 
monotonic increase with respect to 9 as proved by the computer calculation, there is a point p = p n 
where 9 takes n. 

For the trajectory from p v to a distant scattered point, the mcf is written as 



W(r, 9, cf>, E, v, p) = W r -Wg + W^- 2W r (Q) + 2Wg(w). 

The last term has been added to for the continuity of W at the point p — p n to the expression (|4 
The trajectory is given by 



(47) 



— (W r - Wg - 2W r (Q) + 2W e (n)) = — ar g r(a)+V27r, 

d 

<t> = —Wg+ir + (p . 
dp 



(48) 
(49) 
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Fig. 1: A trajectory of an electron of the hydrogen atom in the scattering 'mode' with parameters r] a = 
3.3, v — 2.2, n — 1.1, (solid line) and the classical orbit with 77,, = 3.3, 1 = 2.7, fi = 1.1, (circle). 

The equations of the trajectory prescribed above are more easily seen if the approximate relation 
( ppf ) is used. It shows that coordinate (f> is constant except the point of p v where 6 = it. This indicates 
that the trajectory is almost on a plane like the classical orbit. 

An example of a trajectory of the electron with parameters ~q s = 3.3, v = 2.2, fj, = 1.1,, projected 
on a plane, is shown with the corresponding classical orbit with r/ s — 3.3, 1 — 2.7, /i = 1.1, in Fig. |l|. 

The corresponding classical orbit, as is well known, is very similar to this trajectory except the 
neighborhood of the origin. 

4 Cross section 

Along the line of thought of classical mechanics, the cross section will be detrermined by the mode 
trajectory. 

The differential cross section for the classical orbits of uniform incident beam is given by H 



where 9 is the scattering angle of the electron beam and parameter 77,, (g4|) has been used. This is 
equal to that obtained in quantum mechanics not to mention. The scattering angle is expressed, in 
terms of the angular momentum lh and energy E, as 




(50) 




(51) 



The impact parameter, s, is related to the angular momentum as s = Ihj \j2mE. 



Equation ( (48|) leads to the scattering (or deviated) angle of the electron trajectory, by taking 
p = oo and using the approximate relation (|20|), or d„Wg « 6, 



d d 
= 7r ~ gr,We,p=oc =2— ar g r(a) 



(52) 



In the remote region from the origin, equation ( |37| ) shows that the difference between two positions 
along the trajectory satisfies 



de d d .„ i/ + i/a 



dp <9p dv 



P 



(53) 



Integration gives rise to p9 = z/ + Ya< It is thus obtained that the impact parameter of the trajectory 
is given by 



s = r9 = (u + i/ 2 )h/V2mE. 



(54) 



This indicates that (y+^/^h corresponds to lh, angular momentum in the sense of classical mechanics, 
and v should be greater than — Y2. 
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Fig. 2: Scattering angle as a function of the impact parameter of an incident beam for the classical 
orbit (circle) and for the trajectory (solid line) with ry s = 2. 



An example of the scattering angle of an incident beam as a function of the impact parameter for 
the classical orbit or for the trajectory is shown in Fig. ||. 

The differential cross section for the trajectories of uniform incident beam of the electron is ob- 
tained in the same way as the classical one, by using (p2h, 



<J(0 SC ) 



ds 



de sc 
2 



1 



2mE sin ( 



dp 



(f + Va) 



d 2 

^2 ar S r ( a ) 



d9 s 

-1 



2 sin 6 S 

Parameter v in the right hand side should be expressed in terms of 9 SC through (|52j). 



(55) 
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Fig. 3: Differential cross section as a function of the scattering angle of an incident beam for the 
classical orbit (circle) and for the trajectory (solid line) with ry s = 2. 

An example of the differential cross section of an electron beam for the classical orbit and that for 
the trajectory are shown in Fig. ^. 

The figures indicate that the similarity between the cross sections of the classical orbits and the 
trajectories is, as it were, complete. 

5 Flux 

The usual wave function for the Coulomb scattering is given, by Gordon, as § |ll| 

OO 

ip(r, 9) =J2( 2l + ^ exp(-i argT(l + I + ir h ))Li(r)Pi (cos 9) (56) 
;=o 

where 

(57) 



2' +1 exp(-7r/2r, s ) 

L l( r ) = -ITV I 7 I 7Z \ I U r(P)- 



p\T(l + l + ir) a )\ 



By using (flm and (31), it can be composed of a linear combination of the travelling waves with 



parameters ?] s , v — I and /.t = 

i>(r,6) = Y (21 + 1) (^' "PHH) (^ e -iw r _ Jp re iiw r -2WAQ))) 
^ 7rr(l + / + iris) \ / 

For r large, it is approximately expressed as M, Q 



(58) 



ip(r,9) ~ exp^pcosf? — r/ s logp(l — cos#)] + 



p(l — cos 9) 

x exp [i[p + r] s log/o(l - cos 9) + tt - 2argT(l + irj s )]) . 



(59) 
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The flux (or current) of beam of particles is defined by 



j = J-tp*-X7ip + c.c. (60) 
2m i 



The first term of the right-hand side of the wave function (|59|) gives rise to the incident flux propagating 
along the z axis for r large 

Jin = — z. (61) 
m 

The second term, on the other hand, gives rise to the outgoing flux 

Tik tj 2 

jout = Tf. TiW^- (62) 

The number of particles in the outgoing flux in a solid angle is 

r 2 jp U f r _ fh 2 Tj s \ 2 rfi 
hk/m \me 2 J 4sin 4 y 2 0' { ' 

which gives the differential cross section. 

According to the interpretation envisaged with Gordon's work, the deviated beam should be 
obtained through the scattering of incident particles with various impact parameters. The outgoing 
flux is considered to be an assemblage of particles which, in classical mechanicsl sense, should be 
governed by a plausible causal dynamics. 

In the present treatment, the associated wave with the motion of an electron with rj s , v and fi = 0, 
may be expressed, for r large, as 



= lpv,in + VVut- (64) 

It holds approximately that 

P r ^p- 2 2-^e^\T(l + v + lVs )\ 2 , P g ~ * , (65) 

(2v + 1) smfcJ 

The wave for a beam of particles with various impact marameters should be written as 

rl>(r,0) = ^2aMr,0)u, (66) 



where a v is an arbitrary constant which will indicate the intensity and initial phase of the wave (64). 
The flux is written 

j = 2^ Y^ll ala v ,ii* v -V>ipu> + c.c. = 2£ jV • ( 67 ) 

V v' V v' 

For r large, the incident waves ^Vin's give r i se to the incident flux 

^hk \a u a^\ 7rexp(7T7/ 6 .) |T(1 + v + ir] s )T(l + v' + if] s )\ 



m p 2 sm6 ia a /(2i/+1)(2i/ + 1) 2 »+»'+ 2 
x cos(- argT(l + v + irj a ) + argT(l + v' + irj s ) 
-(0 + l / 2 n)(u - v') + arg(oXO). (68) 
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Similarly, the outgoing waves ?/v,out 's give rise to the scattered flux 

^hk \a v a v i\ 7rexp(7r?7 5 .) |T(1 + v + ir) s )T(l + v' + ir) s )\ 

x cos(argr(l + v + irj s ) — argT(l + v' + ir/ s ) 
+ (9 + y 2 7r)(u - v') - 2tt{u - v') + arg(aX>)) • (69) 

For v ~ i/', 

argT(l + v + ir] s ) ~ argT(l + v' + ii] s ) ~ [y — v')d v argT(l + v + ir/ s ). 

From (§§), 

#out = 7T - (6> scat + 6> in ) = n -2d u argT(l + i/ + ir) s ) - 9 in , 

where 9 m is written as 9 in a preceding section. Thus, the argument of cosine in J^i/'^out with 9 ou ^ 
equals that of cosine in ] UV ',m with 



'in: 



argT(l + v + iri s ) - argT(l + v' + irj s ) + (9 +i/&r)(v - v') 
— 2n(v — v) + arg(a*aj/') 

~ (v - v') [d v argT(l + v + it] s ) + 9 out ~ 3 / 2 tt] + arg(a*cv) 

= -{y - v') [d v argT(l + v + iris) + 9 in - i/a7r] + arg(a*<v). (70) 

It is seen thus that for p ^> 1 

W.in ■ (-f)p 2 sin6» in = j w ', ou t • f /9 2 sin6» ou t. (71) 

If it is assumed that the incident beam \ vv i with v' clearly different from v is incoherent, < 
cos arg(a*a„<) > (average over the initial phases) = 0, the total flux may be expressed as 

= (72) 

v.v' v~v' 

The current \ vv >- with v' ~ v could be discriminated from the one Ja,u' with p! ~ p clearly different 
from v because of the difference of the impact parameter or the angular momentum. Then, expression 
( fn] ) means that the flux density of particles along the trajectory is conserved for r large, or the flux 
density of the incident beam is equal to that of the outgoing beam. Therefore, the cross section for 
the incoherent incident beam can be calculated like the classical mechanical beam as (p5|). 



6 Bound state 

An electron trajectory in a hydrogen atom in a bound state is analyzed. The solution of the form of 
expressions (|l3|) and (jlj) with energy E negative is sought. The mcf in the 9 coordinate space, We, 
is the same in the preceding section, flT7j). 

The traveling wave in the radial coordinate is obtained from Eq. (|^) with replacement of ir) s or 
ip by —77 or p, respectively. By writing u{r, E, v) = rR(r, E,u), a traveling wave is obtained 

r(a)r(6-a) 

u = u v + ru M sinna — — , (76) 

T(b) 

where a = v + 1 — 77, and b = 2v + 2. Functions um and uy are given by 

u M = e~ p p u+1 M(v + l-r 1 ,2v + 2.2p), (74) 
u v = e- p p l,+1 V(v + l-r 1 ,2v + 2,2p), (75) 
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where 



-2mE 
— 9 — 1 



and r\ 



m 

~2~E' 



(76) 



The mcf in the r coordinate space, W r , for the negative energy state is thus given by 



W r (r, E, v) = arctan 



7rr(6-a) M(a,b,2p) 



r(l-o)r(6) V(a,b,2p) 



(77) 



The asymptotic forms of W r at r = or r = oo are summarized in Appendix. 

The mcf for the system is given by Eq. ( |39| ) with expression (fr7|) in place of (|32]). The trajectory 
of an electron in a bound 'mode' (E, v, p) should be given by Eqs. (^) with the same replacement of 
the variables and parameters as mentioned for the scattering state. 

A computer calculation shows that these derivatives dW r /dE, and dW r /dis, are monotonic with 
respect to r. The variation of the derivative dW r /dv between r = and oo cancels out that of the 
derivative dWg/dv between 8 = and tt. It is thus found that the electron moves between r = and 
oo, accompanying with a translation between 8 = and w. 

The trajectory of an electron in a bound state is given as follows. By assuming the origin to 
correspond to the nearest point of classical orbit to the center of the potential and to be the starting 
point of the trajectory, the mcf for the motion from r = to r = oo is given by 



W = W r + Wg + Wfr 



(78) 



where W<j, = p,<f), and W$ or W r is given by Eq. ([D]) or (|77|) ; respectively. The equations of motion 
are, from (0), 



d d 

h dE Wr ^ t = h dE Wr ^ ~ to = ~ <0 ' 

— (W r + We) = —(W r (0) + W 9 (0)) = 0, 
ov av 



(79) 
(80) 
(81) 



where 4>q is the value of <fi at the initial time tg. The mcf for the successive motion from r = oo to 
r = is 



W = -W r -W e + 2W r (oo) + 2W e (n) + W^. 



(82) 



The equations of motion are 



t 



to 



Tl-Ktl d 

l?- h 8E W r 



dv 



{W r + W 6 ) 



0, <j)= — + 



(83) 



by (19), and (A 



in Appendix. 
3 which is just equal to 



where use has been made of d v W r (oo) = —d v Wg(-K) = —n, 

The period of one round trip is found to be —h-nr\l (-E) = 2iu?^Jm/ (-2E) 
that obtained for the classical orbit with the same energy. This suggests that the mcf, (|7?]) with (17), 
is the right one and the mcf could be determined uniquely for the system with the mode parameters. 

A trajectory of an electron in a bound state thus determined is illustrated for the 'mode' with 
parameters rj = 3, v = 2, n = 1, and (f>o = at t = to in Fig. |]. Both of the trajectory or the 
corresponding classical orbit are projected on a plane. It could be recognized that the corresponding 
classical orbit resembles more to the trajectory as the mode parameters rj and v become larger. This 
shows the correspondence principle. 

The requirement that the associated wave ( |l4| ) with the round trip motion of an electron should 
be unique everywhere leads to that each component of the mcf W r , Wg or Wa, must be unique at any 
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Fig. 4: A trajectory of an electron of the hydrogen atom in the bound 'mode' with rj = 3.3, v = 2.2, 
and \i = 1.1 (solid line) and the corresponding classical orbit with I = 2.7 (circles). 



point except multiples of 2ir and parameters 77, v and /1 be integral numbers. It gives rise to the usual 
eigenvalues of the energy for the eigenstates. The wave functions associated with the motion in the r 
and 9 coordinates with rj, v and /1 integral are composed of going and returning waves, Eqs. (|J) and 
(g) with both -2W r (0) + 2W r (oo) and -2W e (0) + 2W e {Tr) being a multiple of 2tt, or Q |l| 

^/W r e lW " - ^/P~ r e- lW - , and y/P~ 9 e lW <> - y/P~ e e- lW <> . (84) 

These are regular for < r < 00 and < 9 < 7r, respectively, and equivalent to the usual stationary 
state wave functions except numerical factors. 



7 Conclusion 

A trajectory of the electron in an attractive Coulomb potential has been shown. The trajectory 
resembles well the orbit of the corresponding state in classical mechanics. It runs also through a 
tunneling region. The trajectory has been derived from the mode characteristic function (mcf) which 
is an extenstion of the Hamilton characteristic function. The period of the round trip motion of the 
electron in the bound state is the same as that of the corresponding classical motion with the same 
mode parameters. The accordance of the characteristics between classical and quantum mechanics 
suggests the validity of the mcf and the dynamics defined. 

It is to be noted that all solutions of the wave equation including waves diverging at singular points 
of the wave equation are necessary to derive the mcf. This suggests the role of all the solutions of the 
wave equation. It has been found that if the wave equation is decomposed into a set of single- variables 
r, 9 and (/>, the resulting differential equations of 2nd order have each mcf uniquely, which should lead 
to a mode trajectory of an electron with some dynamical assumptions. 
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The wave function finite everywhere for any scattering or bound state is made by superposing the 
traveling waves associated with the electron motion, ft is not always continuous for any bound state. 
The requirement that the wave function be continuous or uniqely determined at any point of the 
coordinate is equivalent to that the difference between the mcf's at endpoints be a multiple of ir. It 
leads to that the mode parameters are integral numbers and the eigenvalus of the physical quantities 
become discrete. 

The discussion on the cross section and the flux suggests that the statistical nature of the wave 
theory stems from the fact that the wave consists of waves associated with the particle motion with 
each mode parameter which is causally determined. Not all the characteristics of statistical nature of 
the wave mechanics have been examined but this may be one of the most important characteristics. 

These results indicate the consistent existence of the trajectory in wave mechanics and suggest a 
significance about the relation between a particle motion and a traveling wave function as an extension 
of the de Broglie postulate. 



The author would like to express his sincere thanks to Drs. S. Nakamura, Y. Takano and T. Okabayashi 
for their continual encouragement during the long course of the work. 



A Asymptotic form of W r (a, b, p) 

Asymptotic forms of W r {a,b 1 p) for negative energy states and its derivatives with respect to E and 
v are summarized. It holds from expression ( f77| ) that, by putting F r = tan W r , 



dW r 
~dE 



dW r 
dv 



2E(1 + F?) 



T){ip(b — a) — ip(l — a)} 



_L_ / rjdM pdM \ 1_ fijdV_ pdV \ 
M V da + dp ) + V \~~da~ + ~~dp~ J 



1 + F? 



ip(b - a) + ip{\ - a) - 2-0(6) 



1 fdM 2 9M\ 1 fBV dV 
~M\~da~ + ~db~ J ~ V \~da~ + ~db 



As r tends to 0, it can be written for b > 1 m, [LC 



1 



2a 



M(a, b, 2p) 
V(a,b,2 P ) w T(b-l)[(2p) 1 - b -l] 
By these approximations, it can be obtained that for b > 1 



W r (0) = 0, 



8W r 



(0) = 0, 



8W r 



(0) = 0. 



dE w ' dv 
Asymptotic forms for r large with a and b fixed are [l0| 

M(a,b,2p) « e 2 "(2p) a - b ^|, 
V(a,b,2p) w -e 2p {2p) a - b T{b- a) cos na. 



(A.l) 



(A.2) 

(A.3) 
(A.4) 

(A.5) 

(A.6) 
(A.7) 
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These give rise to 

W r (oo)=-wa, -qe( oc ) = -2E' ^7 (oo) = _7r ' (A ' 8) 
It is found that the variation of dW r /dE between r = and oo is irrj/(—2E). 
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